Microscopic model for spreading of a two-dimensional monolayer 
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We study the behavior of a monolayer, which occupies initially a bounded region on an ideal 
crystalline surface and then evolves in time due to random hopping of the monolayer particles. In the 
case when the initially occupied region is the half-plane X < 0, we determine explicitly, in terms of 
an analytically solvable mean-field-type approximation, the mean displacement X{t) of the monolayer 
edge. We find that X{t) ~ A^/D^i, in which law Dq denotes the bare diffusion coefficient and the 
prefactor A is a function of the temperature and of the particle-particle interactions parameters. 
We show that A can be greater, equal or less than zero, and specify the critical parameter which 
distinguishes between the regimes of spreading (A > 0), partial wetting (A = 0) and dewetting 
(A < 0). 



1 Introduction 

It has been well appreciated since the pioneering work by Hardy fl] that a liquid droplet spreading on a 
solid, although it undergoes no visible change in shape, emits a very thin invisible film - the precursor, 
which advances at a seemingly faster rate than the nominal contact line. Hardy was able to detect 
its presence by observing a significant change in the value of the static friction of the surface. Stating 
that he was unable to conceive of a mechanism by which the film can be emitted and spread further 
along the solid substrate, Hardy proposed that spreading of the film occurs by a process involving a 
continual condensation of vapor. Bangham and Saweris Q have demonstrated, however, that such a 
film shows up even in the absence of any vapor fraction, suggesting thus that the physical mechanism 
giving rise to the precursor film can be different from the evaporation/condensation scheme. Later, 
using ellipsometric and interferometric techniques, Bascom et al. @ investigated the spreading of 
the precursor film from a more quantitative point of view. They examined the behavior of various 
non-polar liquids on clean metal surfaces in the presence of both saturated and unsaturated air and 
concluded that the precursor film is always present; making the air saturated or unsaturated with 
vapor, roughening the surface and purifying the liquids does not eliminate the film, but only affects the 
speed at which it spreads over the solid. The thickness of the film was also found to depend essentially 
on the liquid/solid system under study; it can be as small as molecular size (several angstroms) or can 
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amount to hundreds of angstroms (see [|j for a review) . Particularly, the droplet of squalane spreading 
on stainless steel exhibited a precursor with a thickness of approximately 20 angstroms. 

Spreading rates and dynamical shapes of advancing precursor films have been studied thoroughly 
for many years, both experimentally and theoretically (see [|, || fj], ||] and references therein). 
These studies have resulted in a rather good understanding of the problem. It was realized that the 
precursor film appears not because of the condensation of the vapour fraction, which process may, of 
course, exist but plays a minor role. Rather, such a film extracts from the droplet and advances along 
the solid surface due to the presence of attractive interactions between the fluid molecules and the solid 
atoms [||, H, pL o, m. Later, it even became possible to elaborate consistent hydrodynamic theory of 
spreading of thin precursor films, based on the celebrated lubrication approximation of fluid mechanics 
[@> ||) |6|) Hi- in particular, this considerable theoretical advancement allowed to resolve an old-standing 
enigma concerning the dynamics of the nominal contact line: it was well documented experimentally 
(see, e.g. || and references therein) that in the complete spreading regime the radius R of the 
nominal contact line shows a slow growth with time, R ~ i 1 / 10 , where the exponent 1/10 is universal, 
i.e. dependent only on the geometry and independent of the precise nature of the liquid/solid system. 
As a matter of fact, in this law even the prefactor appears to be insensitive to the spreading power 
S, which is the characteristic of a given liquid/solid system and equals the difference of the surface 
tensions of the solid-vapour, solid-liquid and liquid-vapour interfaces respectively. On the other hand, 
conventional analysis, which considers the gradient of the liquid/solid free energy as the driving force 
of spreading, predicts much faster growth which is, moreover, dependent on the spreading power. The 
hydrodynamic picture developed in [||, ||] has found the solution of this controversy, showing that 
the spreading power S is totally dissipated into the precursor film, such that the dynamics of the 
nominal contact line is actually not affected by S. Besides, it has been recognized that contrary to 
the general belief, the final state in complete (S > 0) spreading of a liquid droplet is not necessarily 
a monolayer covering the solid surface: it is only the macroscopic part of the droplet which spreads 
completely; the precursor film, however, can cease to spread further after the macroscopic part of the 
droplet is exhausted and form a stable, nearly planar droplet-like structure - the so-called "pancake" 
H |{| [l(| 11, |l~2l| . The thickness of the pancake is fixed by a competition between S and long-range 



attractive interactions with the substrate, which tend to thicken the film. For pure van der Waals 
interactions the thickness e p (S) of the pancake is e p (S) ~ a(3j/2S) 1 ^ 2 , where a is the molecular size 
and 7 is the macroscopic surface tension of the liquid. In general, 7/5 ~ 1 and e p {S) ~ a, but if S is 
small, e p (S) can be relatively large. Therefore, within the framework of the hydrodynamical approach 
it is the pancake thickness e p (S) which sets the upper bound on the thickness of the precursor film. 
Details concerning different "pancake" structures and their experimental observations can be found in 
a recent review [Ol. 



1.1 Experimental studies of spreading of molecularly thin films. 

Hydrodynamic approach, however, presumes a certain lower cut-off length comparable to the molecular 
size, below which it is not justified. With the advent of modern experimental techniques, capable of 
studying behavior of precursor films whose thickness is in the molecular range, it has become clear that 
the developed theoretical concepts concerning the dynamics and equilibrium properties of thin wetting 
films can only aplly to sufficiently thick films. In fact, for molecular films significant departures from 
the hydrodynamic picture have been observed experimentally [14, [l6|, |l7], |l§|| . 

Extensive ellipsometric measurements [14, 15, pi], |l7], [l^] have elucidated several remarkable fea- 
tures, which seem to be quite generic for spreading molecularly thin films: 
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(i) Experimental measurements carried out on different substrates and with different types of sim- 
ple liquids, polymeric and surfactant melts, showed unambiguously that the radius of the precursor 
film grows with time t as i 1 / 2 . The exponent 1/2 appears to be completely independent of the nature 
of the species involved. The latter affects the spreading rate only through the prefactor in the yt-law. 
Essentially the same behavior has been discovered in the capillary rise geometry, in which a vertical 
solid wall is put into contact with a bath of liquid. The length of the molecularly thin film extracting 
from meniscus and climbing upwards along the wall was also shown to grow in proportion to the 
square-root of time [17]. This shows that the y^-law is seemingly independent of the precise curvature 
of the film's edge, which is a circular line in case of the sessile drops and planar in the capillary rise 
geometries. 



(ii) It was found that the particle density along the film is not necessarily constant. In several 
experimental situations, particularly, for the droplets of squalane, the radial density was seen to vary 
strongly with the distance from the nominal contact line and the variation was progressively more 
pronounced for larger spreading rates. 



(iii) For intermediate-energy substrates a fascinating transient regime of " terraced spreading" was 
discovered [14, [15]], in which several molecularly thin precursor films extract from the macroscopic 
droplet and spread on top of one another each one growing as y/i. Such a regime may last within 
a considerably long period of time, until the layers on the top run out and one eventually ends up 
with a bounded monolayer on the solid surface. Depending on the physical conditions, the monolayer 
may then either continue to spread with the radius growing in proportion to (but with a different 
prefactor, of course), or may remain in form of a wetted spot [17]. The remarkable "terraced wetting" 
effect has been also seen in the capillary rise geometries. Besides, recent numerical, Molecular Dynam- 
ics simulations of liquid drop spreading were able to reproduce both the "terraced" wetting spreading 
regime and the y/i-law (see p9|, 30, 31] and references therein). 



1.2 Theoretical studies of spreading of molecularly thin films. 

Theoretical analysis of the physical mechanisms underlying the seemingly universal y/i-law and the 
"terraced wetting" phenomenon followed three different lines of thought: 

In an analytical description of the "terraced wetting" phenomenon has been proposed, in which 
the liquid drop on a solid surface was considered as a completely layered structure, each layer being 
a two-dimensional (2D) incompressible fluid of molecular thickness but with a macroscopic radial 
extension. The interaction energy of a molecule in the nth layer with the solid substrate was taken 
in the general form as a negative, decreasing function of the distance from the substrate. Next, it 
was supposed that spreading proceeds by filling the successive layers by the molecules from the above 
layers, which process is favored by the attractive interactions with the solid. In each layer there is a 
horizontal, radial particle current and vertical permeation fluxes, one from the upper layer and one 
towards the lower layer, which appear to be located in a very thin "permeation ribbon" just at the 
droplet surface - the core of the droplet seems to be a "stagnant" liquid with respect to the vertical 
mass transfer. In such an approach, de Gennes and Cazabat [|E| have found that whenever the distinct 
layers grow at a comparable rate, they grow in proportion to t 1 / 2 . In case when the precursor film 
grows at a much faster rate than all other layers, the model predicts that the precursor grows in 
proportion to (t//n(t)) 1 / 2 . Apart from the theoretical prediction of the y/i-l&w, the model allows to 
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draw certain conclusions concerning the appearence of the "terraced wetting" regime. To some extent, 
the underlying assumptions and predictions of this model concerning the intermediate time behavior 



were confirmed by Molecular Dynamics simulations [20] 



An alternative description within the framework of non-equilibrium statistical mechanics has been 



elaborated in [21, 22, 23]. Here, an interfacial model for the non-volatile fluid edge has been developed 
and analysed in terms of Langevin dynamics for the displacements of horizontal solid-on-solid (SOS) 
strings {hj} at increasing heights j = 1,...,L from the substrate, which strings have essentially the 



same meaning as the layers in the de Gennes-Cazabat model [19]; hj can be thought of as the radius 



of the jth layer. Anticipating the discussion of results obtained in the present paper, we will describe 



this promising theoretical approach giving some more details. In the model developed in [^, 22, 23] 
the energy U({hj}) of a given configuration {hj} was described by 

L 

U(h ,hi,...,h L ) = ^P{hj - hj-i) - no h , (1) 
i=i 

where ho is the length of the "precursor" film, fio is the wall tension and the function P{hj — hj^i) 
determines the contribution due to the surface energy. Explicitly, it was taken as 



P(hj - hj-!) = J yl + (hj - Vi) 2 , (2) 

where the parameter J is related to the surface tension; thus, the spreading power S for such a model 
is given by S = /io — J- In the limit hj — hj-\ <C 1, the function in Eq.(2) reduces to the standard 
Gaussian form, i.e., 

P(hj-hj„ x ) « Jihj-hj^) 2 , (3) 

while for \hj — hj—i\ 3> 1 it obeys 

P(hj-hj-i) ps J \hj - hj-i\ (4) 
In other words, Eq.(2), which bears a certain relationship with the Lifschitz equation describing the 



time evolution of an interface due to the effects of the surface tension |2J], supposes the following, quite 
realistic behavior of the interfacial energy: for small distortions of the interface the surface energy is 
Gaussian, while the cost of a large distortion is linear with the distortion size. Dynamics of the strings 
{hj} is then described by the set of L coupled Langevin equations 

dh k dU({h 3 }) 
t~dt = - —dh^ + f{hk ' t] > (5) 

where £ denotes the friction coefficient, which is supposed to be the same for all layers, and f(hk',t) is 
Gaussian, delta-correlated noise. The model in Eqs.(l),(2) and (5) allows for an analytical, although 
rather complicated solution, which shows an extraction of the precursor film and "terraced" forms 
of the dynamical thickness profiles. It predicts that for fiQ > J, (S > 0), and for sufficiently short 
precursors, (such that the dominant contribution to the surface energy is given by Eq.(3)), the length 
of the film increases with time as 

ho(t) - sj [fio - J) t, (6) 

which resembles the experimentally observed behavior. For sufficiently large precursor, for which the 
energy increases linearly with the length, Eq.(4), the layer on top of substrate is found to show a faster 
growth 

ho(t) ~ (no~J)t (7) 



Finally, it was found that exactly at the wetting transition point fig = J, the precursor film advances 
in proportion to 



ho(t) ~ y/t ln(t) (8) 

Consequently, this model predicts that at very large times the advancing precursor film attains a 
constant velocity; the \ft-\am is thus found only as a transient stage. Besides, dynamics of the layers 
at large distances from the substrate and, respectively, of the macroscopical dynamical contact angle 
disagree with experimental data. Apparently, this inconsistency with experiments can be traced back 
to the fact that focusing on the evolution of the interface only, the model neglects dynamics in the 
liquid phase and thus discards the energy losses due to viscous flow pattern, generated in the spreading 



droplet. In [21, g^, |2J| a viscous- type dissipation is assumed with constant friction coefficient, which 
represents rather generic and oftenly loose assumption used in the descriptions of the phase-separating 
boundary dynamics in terms of the time-dependent Landau-Ginzburg-type model. In fact, the model 
underestimates the dissipation in each layer hj. Account of the energy losses in viscous flows in the 
core region of the droplet and of the dissipation in the vicinity of the solid substrate actually results 
in the overall viscous-type dissipation, which is not surprising for the system with many degrees of 
freedom. The friction coefficient, however, turns out to be weakly dependent on the height above 
the substrate and, what is essentially more important, appears to be an increasing function of hj. 
Therefore, these dissipation chanels should be certainly taken into account within the framework of 



the powerful theoretical approach proposed in [g^, 22, 23|, which may result in a consistent dynamical 



theory of partial and complete wetting valid for all scales. Such improvements are currently under 



investigation [25| 



Lastly, a microscopic dynamical model for spreading molecularly thin films has been devised in 
[p6| , pTj ]. Here the film was considered as a two-dimensional hard-sphere fluid with particle-exchange 
dynamics. Attractive interactions between the particles in the precursor film were not ostensibly 
included into the model, but introduced in a mean-field-type way - it was supposed that the film is 
enclosed by the SOS-model interface, in which the parameter J was treated as some (not specified in 



[|26| , 27]) function of the amplitude of the particle-particle attractions. The film was assumed to be 
connected to a reservoir of infinite capacity - the macroscopic drop. The rate at which the reservoir 
may add particles into the film was related to the local particle density in the film in the vicinity 
of the nominal contact line and to the strength of the van der Waals attractive interactions between 
the fluid particles and solid atoms within the framework of the standard Langmuir adsorption theory. 



Contrary to [19| and to the hydrodynamic picture of ||], the model in |26|, 27] emphasized the issues of 
compressibility and molecular diffusion at the expense of the hydrodynamic flows; it was assumed that 
the reservoir and the film are in equilibrium with each other, so that there is no flow of particles from 
the reservoir which pushes particles to move along the substrate away of the droplet. In this approach, 
the y/i-law for growth of the film was first analytically obtained for the capillary rise geometries and it 
was actually found that the density in the film does varies strongly with the distance from the reservoir. 
This agrees, at least qualitatively, with experimental observations (see the Introduction, 1.1, (ii)), but 
no direct comparison was made, as yet. Clearly, the factor which makes such a comparison quite 
awkward is that the spreading rate is expressed through the parameter J, which is supposed to be 
some known parameter. As it will be made clear below, this parameter depends on the particle-particle 
attractions and moreover, on the density distribution in the precursor film. Further on, the critical 
conditions under which spreading of the precursor film may take place were established in p6l [27]]. 
It was also suggested that the physical mechanism underlying the \/t-l&w stems from diffusive-type 
transport of vacancies from the edge of the advancing film to the macroscopic liquid edge, where they 
perturb the equilibrium between the macroscopic drop and the film and get filled with fluid particles 



from the macroscopic liquid drop. In |2q] this picture was extended to the case of sessile drops and 



it was shown analytically that the curvature corrections result only in a weak slowing down of the 
precursor spreading; the film radius grows in this case as {t/ln(t)) l l 2 , which prediction agrees with 



To close this introductory part of our paper we mention several analytical studies of the process, 
which can be thought of as the reverse counterpart of wetting, - dewetting of microscopically thin 
liquid films from solid substrates. Recently Ausserre et al. investigated analytically dewetting 
of a monolayer, which was assumed to proceed by nucleation of holes (bare regions) and creation of 
"towers" - two-layer regions. Considering the monolayer as an incompressible 2D liquid, it was shown 
that the hole radius R (or the radius of a "tower") grows with time in proportion to {t/ln(t)) 1 / 2 . 
This means that in the monolayer regime the dewetting process proceeds essentially slower than in 
the case of mesoscopically or macroscopically thick films, for which the behavior R ~ t is generic 



3], gj, 35, 36 1 . Another interesting example of a (forced) dewetting of a monolayer was discussed in 



and concerned with the squeezing of a molecularly thin liquid film out of a narrow gap between 



two immobile solids. In [37| the mechanism responsible for squeezing was attributed to the process of 
spontaneous opening of holes in liquid layers; the holes are subsequently get filled by deforming solid 
material exerting pressure on the hole boundaries. Viewing solids as isotropic, structureless elastic 
media and the liquid phase in a lubricated contact between two solids as a sequence of layers, each 



layer being an incompressible 2D liquid, Persson and Tosatti [37] were able to estimate the critical 
radius of the hole, necessary to initiate further squeezing, and to define the rate of the removal process 
after the nucleation of a critical hole has occured. It was shown that the radius of the hole, which 
exceeds initially the critical value, grows in time again in proportion to (t/ln^t)) 1 ^ 2 . 



2 The objectives and a brief outline of the paper. 

In this paper we study analytically, in terms of a stochastic microscopic model, the behavior of a 
liquid monolayer in a situation, in which a monolayer occupies initially only some part of solid surface 
- the half-plane X < 0, Fig.l, and then is allowed to evolve in time due to the thermally activated 
random motion of the monolayer particles. Here we aim to calculate the time dependence of the 
mean displacement of the monolayer edge, defined as the position of the rightmost monolayer particles 
(Fig.l), to determine the prefactor in this dependence in terms of the interaction parameters and the 
edge tension of the monolayer. 

We note that the model to be considered here clearly shares common features with many dewetting 
and wetting experimental situations and models, which were described in the Introduction. A mono- 
layer in such a non-equilibrium configuration appears, for instance, at the late stages of spreading, 
when the liquid drop feeding the precursor film gets exhausted or in the situation when the monolayer 
on the solid surface is perturbed by a sudden removal of some amount of particles or by nucleation of 
a dewetted region - a circular hole or a patch. Such a model applies, after some minor modifications, 
to the dynamics of ultrathin liquid columns in nanopores or in narrow slits between solid surfaces 
(Fig. 2). It can also serve as a microscopic description of the process of Ostwald ripening of voids, 
spinodal decomposition or island formation in two-dimensional adlayers. We note also that in view 
of the above-mentioned results concerning spreading dynamics, we expect that the precise geometry 
of the two-phase region is not very important; the difference between the case when the front of the 
monolayer is planar, as we consider here, and the case when it is a circular closed line, what should 
be for circular dewetted holes, can be only in the appearence of logarithmic in time corrections to the 
\/t-law, important at times when the displacement of the edge becomes comparable to the radius of 
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the edge curvature. An important point is that both the region occupied by the monolayer and the 
initially dewetted region should be macroscopically large. 

As opposed to [19], [32] and |J7J, we will regard the monolayer as an essentially discrete, molecular 



liquid composed of interacting particles moving randomly on an ideal crystalline surface. Particles 
migration is assumed to be activated by the solid atoms vibrations and will be described using the 
standard Kawasaki picture for particles exchange dynamics under long-range particle-particle interac- 



tions. The picture we make use thus follows closely the model elaborated in [26, 27], being different 
from the latter in two important aspects: first, the long-range attractive interactions between the 
liquid particles are here explicitly included into the dynamics, and, consequently, our results will be 
expressed in terms of the interaction parameters. Second, the reservoir of particles is absent, which 
allows us to study within a unified approach dynamics of both spreading and dewetting processes. We 
also hasten to remark that with regard to other dynamical wetting theories, this model is related to the 



Molecular Kinetic theory of wetting dynamics, proposed and developed by Blake et al. [38, 38 , 4C]. In 
this theory, which emphasizes the dissipation in the vicinity of the nominal contact line at expense of 
the dissipation due to viscous flows in the "bulk" droplet, the analysis of dynamics of spreading liquid 
droplet was reduced to a mean-field-type consideration of the forced, thermally activated motion of 
fluid particles which appear directly at the droplet edge. In our case, however, the driving force is not 
assumed a priori, but is found consistently as the result of the cooperative behavior, associated with 
the interplay between the long-range attractive particle-particle interactions and repulsion at shorter 
scales. As well, we deal here with simultaneous random motion of all particles in the film, not reducing 
the problem to consider the dynamics of only particles at the edge. 

Next, we will make several simplifications, compared to real physical systems. In what follows 
we will assume that creation of "towers" and particles evaporation in the direction normal to the 
solid surface are completely suppressed by the liquid-solid attraction, and hence we will constrain 
our consideration to the system which always remains in a two-dimensional world. Shortcomings of 
this picture will be discussed below. These simplifications will permit us to focus exclusively on the 
dynamical processes which take place in two-dimensions and thus to single out the behavior which 
stems from the interplay between the compressibility and the intermolecular interactions. We note also 
that these mechanisms are entirely complementary to those discussed in [^, 32, 37] and, consequently, 



understanding of their impact on the monolayer evolution is necessary for a complete picture of the 
phenomenon. We remark also that such an assumption can be clearly relaxed for liquids in confined 
geometries, e.g. in nanopores or in narrow slits between solids, where the geometry itself rules out an 
appearance of the vapour phase and thickening of liquid films (Fig. 2). Here we will focus, however, 
solely on the situation with a monolayer on top of open solid surface; relevant cases, as depicted in 
Fig.2, will be discussed elsewhere. 

The paper is outlined as follows: In section 3 we describe the model and write down basic equations. 
In section 4 we discuss an approximate approach to the solution of dynamical equations. Section 5 
presents the results. Finally, in section 6 we conclude with a brief summary of our results and 
discussion. 



3 The model and basic equations. 

We proceed further with more precise definitions related to the model to be studied here (see also 
[p6| , |27f[ for a detailed discussion). Particles of the monolayer experience two types of interactions; 
interactions with the solid atoms (SP) and mutual interactions with each other (PP). The SP interac- 
tions are characterized by a repulsion at short scales and a weak attraction at longer distances. The 
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SP repulsion keeps the monolayer particles at some short distance apart from the surface, while the 
attractive part of the SP potential hinders particles desorption. Following p6l |27| we assume here 



that the SP interactions correspond to the limit of the intermediate localized adsorption JO], 42]: the 
monolayer particles are neither completely fixed in the potential wells created by the SP interactions 
(Fig.l), nor completely mobile. Potential wells are very deep with respect to desorption (desorption 
barrier 3> kT) so that only a monolayer can exist, but have a much lower energy barrier V\ against 
the lateral movement across the surface, Ud 3> V\ > kT. In this regime an adsorbed particle spends 
a considerable part of its time at the bottom of a potential well and jumps sometimes, due to the 
thermal activation, from one potential minimum to another. Thus, on a macroscopic time scale the 
particles do not possess any velocity. 

We note that such a type of random motion is essentially different of the standard hydrodynamic 
picture of particles random motion in the two-dimensional "bulk" liquid phase, e.g., in free-standing 
liquid films, in which case there is a velocity distribution and spatially random motion results from 
the mutual particle-particle interactions; in this case the dynamics can be only approximately con- 
sidered as an activated hopping of particles, confined to some effective cells by the potential field of 
their neighbors, between a lattice-like structure of such cells (see, e.g. [43, 44]). In contrast to the 



dynamical model to be studied here, standard two-dimensional hydrodynamics pressumes that the 
particles do not interact with the underlying solid. In realistic systems, of course, both the particle- 
particle scattering and scattering by the potential wells due to the interactions with the host solid, as 
well as the corresponding dissipation, are crucially important [45 , 46]; the latter, particularly, remove 



the infrared divergencies in the dynamic density correlation functions and thus make the transport 
coefficients finite [47, |4q ]. Complete dynamical description of particles migration on the solid surface 



can be approached apparently along the lines proposed in [47, |4£[] or, on a microscopic level, in terms 
of the cellular automaton-type description of |4jJ; here we will be thus concerned only with a certain 
approximate model of particles dynamics, appropriate for situations in which the particle-particle in- 
teractions are essentially weaker than the particle-solid interactions. We note that such an assumption 
actually makes sense since the latter are usually at least ten times greater than the PP interactions, 
and therefore should not be appreciably affected by the lateral interactions of adsorbed particles. 

Turning next to the particle-particle interactions, we suppose that these are additive and central, 
i.e. the interaction potential U(fj,ri) depends only on the distance r = \ fj — ri\ of separation of the 
jth and ith particles, U{fj,ri) = U(r). Particles are assumed spherical so that no orientation effects 
enter and we take that the potential energy between a pair of adsorbed molecules is given by 

, + oo for r < a 

U(r) = (9) 
-U (T)(a/r) 6 for r > a, 

i.e. we use the "hard-sphere" core and the usual r -6 attractive term for large r; the minimum occurs 
at r = a for which U(r = a) = —Uq(T). The argument (T) in the parameter Uq(T), Uq(T) > 0, 
signifies that in general case, this property can be dependent on the temperature. Particularly, for 
the Keesom-van der Waals interactions one has Uq{T) ~ 1/T. For the London- van der Waals particle- 
particle interactions Uq{T) does not vary with the temperature. As we have already mentioned, we 
will suppose in what follows that the amplitude of the particle-particle attraction Uq(T) is less than 
the barrier for the lateral motion, i.e. Vi, such that the particle-particle interactions do not perturb 
significantly the array of potential wells created by the particle-solid interactions. 

Now, we specify the particle dynamics more precisely (see also [^, ^7|). Under the physical 
conditions as described above, we can regard the particles dynamics on the solid surface as an activated 
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hopping between the local minima of an array of potential wells, created due to the SP interactions 
PH , f42[ . Thus particles' migration on the surface proceeds by rare events of hopping from one well 
to another in its neighborhood. The hopping events are separated by the time interval r, which is 
the time a given particle typically spends in each well vibrating around its minimum; r is related 
to the temperature T = the barrier for lateral motion V\ and the frequency of solid atoms 

vibrations u) through the Arrhenius formula. We thus may estimate the diffusion coefficient for such 
a motion (which will be exactly the diffusion coefficient of an isolated particle on the solid surface) as 
Dq ~ l 2 /zt, where z is the coordination number of the lattice of wells and I is the interwell spacing. 
In what follows we will suppose that I ~ a, i.e. that the radius of the particle-particle hard-core 
and the interwell spacing are approximately the same. Now, diffusion coefficient Dq will be the only 
pertinent parameter describing the evolution of the local density in the monolayer in absence of the 
PP interactions. When the latter are present, as we actually suppose here, dynamics of any given 
particle is fairly more complicated and is coupled to the instantaneous configuration of the monolayer 



particles (see, e.g. [50, 51, [5^, for discussion). That is, for any particle, releasing from the well 
with radius-vector r, not all hopping directions are equally probable and the particle has a tendency 
to follow the local gradient of the energetic surface U(r;t), created by the mutual PP attractions. On 
the other hand, hard-core repulsion imposes sterical constraints preventing particles crossing and thus 
the double occupancy of any potential well. More specifically, we will account for the PP interactions 
as follows: we will suppose that releasing at time moment t from the well with radius-vector r any 
given particle first "choses" the direction of jump with the (position- and time-dependent) probability 

p(r\?) = Z- 1 exp (J [U(r;t) -U(P;t)]) , (10) 

where r' is the radius- vector of one of z wells neighboring to the well at position r, Z is the normalization 
factor, defined as 

'P [~T T / — * 



Vexp 



U(r;t)-U(r';t) , (11) 



in which the sum runs over all wells neighboring to the well as position r, and the PP energy landscape 
is determined by 



U(F;i) = -U (T)a»J2^pJt (12) 



r 



In the latter equation the summation with respect to r" extends over the entire surface, excluding 
r" = r, and r](r",t) is the time-dependent occupation variable of the well at position r" at time t; 
rj(r",t) = 1 if the well is occupied by a monolayer particle and rj(r",t) = if it is empty. 

Finally, hard-core part of the interaction potential in Eq.(9) will be taken into account in the 
following way: we suppose that when the jump direction is chosen, the particle attempts to jump into 
the target well. We stipulate, however, that the jump can be only then fulfilled, when at this time 
moment the target well is empty; otherwise, the particle attempting this hop is repelled back to its 
position. 

In such a picture of particles dynamics and interactions, which represents, in fact, the standard 
formulation of a hard-core lattice gas dynamics under long-range particle interactions, the evolution 
of the local occupation variable r](f; t) can be described by an appropriate probabilistic generator 



L{r](r; t)} (see, e.g. |51, |52|]). Here we will not go into the details of rigorous probabilistic formulations, 
and will proceed by making a simplifying physical assumption that the realization average of the 
product of the local occupation variables of different wells factorize into the product of their average 
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values, which corresponds to the assumption of local equilibrium. It was shown recently in [54, 55] 
that such an assumption provides an adequate description of particles dynamics in hard-core lattice 
gases and we thus expect that it will be also a fair approximation for the system under study. The 
assumption of the local equilibrium allows us to describe the system evolution in terms of local densities 
p(f; t), p(f; t) = r/(r; t), which define the probability of having at time t a particle in the well at position 
f. Consequently, instead of the probabilistic equations describing evolution of rj(r;t), we will have to 
consider a deterministic integro-differential equation describing evolution of the the local densities 
p(r;t). In doing so, we find then that the dynamics of p(r;t) is governed by the following continuous- 
time equation 

dp(f ; t) 



dt 



p(r;t) (l - p(?;tj) 



+ 



+ (1 - p(f;t)) Vp(r'|r) p(r';t), (13) 



where the realization-average transition probabilities are given by 



p(f\r') = Z 1 exp 



pU (T)o* 



p(r"-t) ^ p(r";t) 



(14) 



Eq.(13) has a simple physical meaning - it describes the balance between the departures of a particle 
from the well at position r to any of the neighboring wells and the arrivals of particles from the 
neighboring wells to the well at position r. Particularly, the first term on the right-hand-side of Eq.(13) 
describes all possible events in which a particle, occupying at time t the well at r (the factor p(r;t)) 
may jump, at a rate p(r\r') prescribed by the corresponding change in the energy of the monolayer, to 
any of vacant (the factor (1 — p(r';t)) adjacent wells. In a similar fashion, the second term describes 
the corresponding (positive) contribution due to arrivals of particles from adjacent wells to the well 
at position r. 

Eq.(13) has to be solved subject to the initial condition 

jo for X > 

p(r;0) = { (15) 

I p for X < 0, 

where p denotes the initial mean coverage (number of occupied wells per total number of wells in a 
unit area) of the half-plane X < 0. Eqs.(13) and (15) represent the mathematical formulation of the 
problem under study and allow for the computation of the monolayer edge time evolution. 



4 Approximations. 

One possible approach is to seek for an approximate solution of Eqs.(13) and (15), turning to the 
continuous-space limit and expanding the local densities into the Taylor series up to the second order 
in powers of a and the exponentials in Eq.(14) up to the first order in the gradient terms. In doing 
so, we obtain from our Eq.(13) the following continuous-space Fokker-Planck-type equation with non- 
local, configuration-dependent potential term 

= D [Ap(f;t) - P U (T) V {p(f;t) x 

x (1 - p(r;t)) [ dr' p{r';t)V 1 }] (16) 

J \f — T \ 
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Equation (16), which was rigorously derived by Giacomin and Lebowitz [52| for related lattice-gas 
model with Kac potentials, allows for an analytic, although rather complicated analysis (see [52] for 
more details). 

Here we will pursue, however, more simple approach of |2(| ^7j, which allows for a quite straight- 
forward computation of the mean displacement of the monolayer edge and of the edge tension directly 



from Eq.(13). Following [26, |27j] we assume that for the long-range, but rapidly vanishing interaction 
potentials as defined in Eq.(9), a hop of any monolayer particle being in the "bulk" monolayer does not 
change the energy in Eq.(12). It means, in turn, that for such a particle all hopping direction appear 
to be equally probable and the hopping events are constrained by the hard-core interactions only. This 
is certainly not so for the particles being directly at the edge of the monolayer; these always have "free 
space" in front of them and the monolayer particles behind, which will result in effective attraction of 
the edge particles to the "bulk" monolayer and, in consequence, in asymmetric hopping probabilities: 
the edge particles will attempt preferentially to hop towards the "bulk" monolayer. In reality, all 
particles of the monolayer will experience a weak attraction in the negative X-direction. Being zero 
for X = — oo, the effective attractive force will grow and reach its maximal value for X = X(t), i.e. 
for the edge particles. Simplifying the actual picture to some extent, we will suppose here that this 
"restoring" force is present only for the particles which are directly at the edge (see also p6| , p7|). This 
resembles, in a way, the model in [21, p2[ 23|, which was concerned with the dynamics of the edge due 



to the edge tension only. In contrast to this model, however, the approach of |26|, 27] does not neglects 
the presence of the "bulk" monolayer phase; as we will see below the hard-core interactions between 
the particles in the "bulk" monolayer are crucially important and, particularly, are responsible for the 
\/t-behavior in place of the linear in time dependence predicted by Eq.(7). 

In a more precise way, the main assumptions of the approach in |2q, ETfl can be formulated as follows: 



(a) Assume that at any time moment the monolayer is homogeneous in the direction normal to 
the X-axis. Consequently, the energy U (r ; t) stays constant for any hop which does not change the 
particle position along the X-axis and the probability for such hops to take place is site-independent 
and is equal to 1/z. This implies, in turn, that the edge of the monolayer is sharp (straight line) 
and allows to reduce the problem to the effectively one-dimensional model, in which the presence of 
the second direction is accounted for only through the renormalized diffusion coefficient and actual 
two-dimensional tension of the monolayer boundary. 

(b) In the general case the monolayer edge will change its position along the X-axis, i.e. the 
monolayer will either contract or dilate in the X-direction, and thus the density will be dependent on 
the X-coordinate. One may suppose, however, that the density distribution p(X; t) will be a slowly 
varying, at a microscopic scale, function of X, such that for weak long-range potentials in Eq.(9) and 
for X which are strictly less than the instantaneous position of the monolayer edge X(t) the difference 
U(X + a; t) - U(X - a; t) < 1/(3. 



(c) Hopping probabilities of particles being directly at the edge obey Eq.(14) and thus depend 
implicitly on both X(t) and p(X;t). On the other hand, one may expect that after some transient 
period of time (which will be not studied here), the probability of making a jump away of the edge, 
i.e. p(X(t)\X(t) + a), and the probability of making a jump towards the "bulk" monolayer, i.e. 
p(X(t)\X(t) — a), approach some limiting values p and q, which do not depend on X(t) and t. We 
note that this expectation is actually consistent with the solid-on-solid-model description of the liquid- 
vapour interface [21|. We will show below that this is actually the case and stems from the stabilization 
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of the density profile in the vicinity of the moving edge. 

Finally, p and q are to be determined in a self-consistent way. To do this, we will first solve the 
model described by (a) - (c) supposing that p and q are known, fixed parameters, and calculate the 
time evolution of X(t) and p(X,t). Then, substituting the latter into Eq.(14), we will obtain the 
closed- form transcendental equation, which determines the dependence of p and q on Uq(T), p and ft. 



5 Results. 

Let us consider now the approximate picture of the monolayer evolution, described by (a) - (c), 
supposing first that p and q are some given parameters. We notice that such a picture has an evident 
interpretation within the framework of a lattice gas dynamics (Fig. 3). It defines, namely, the evolution 
of a symmetric hard-core lattice gas, which is initially placed, at mean density p, at the sites — oo < 
X < of a one-dimensional infinite lattice of regular spacing a. The particles are allowed to perform 
random hopping motion between the nearest-neighboring sites; the hopping motion is constrained by 
hard-core interactions. All particles, excluding the rightmost particle of the gas phase, have symmetric 
hopping probabilities, i.e. for them an attempt to hop to the right and an attempt to hop to the left 
occur with equal probability m = 1/zt. In contrast, for the rightmost particle these probabilities are 
asymmetric; they are equal to p and q for jumps away of and towards the gas phase respectively. 

Dynamics of the rightmost particle in such an asymmetric, (with respect to the density distribution 
and the hopping probabilities of the rightmost particle), lattice-gas model has been analysed recently 
in ]55| ]. It was shown that for arbitrary values of the ratio p = p/q, (0 < p < 1), the mean displacement 
of the rightmost particle follows 

X(t) = A VA^, (17) 

in which equation Do stands for the diffusion coefficient, Dq = a 1 /zt, and the parameter A is deter- 
mined implicitly as the solution of the transcendental equation 

^ exp(A 2 /4) [1 + HA/2)] = 11 = (18) 
2 1 — fi 

where $(x) denotes the probability integral. Besides, it was shown in [55] that at sufficiently large 
times the density distribution past the rightmost particle obeys 

P(A;i) = l+ P I{A) I 1 + A2 J" dz exp(-^-(z 2 -2z))}, (19) 

in which 

1(A) = ^Aexp(A 2 /2) [1 + ^A/y/2)], (20) 

and 9 is the scaled variable, 6 = A /Ay/ Dot, where A stands for the relative distance from the rightmost 
particle, A = X(t) - X. In the limit A < y/Dtf/A, Eq.(19) reduces to 

p(X;t) » (1 - M )[l + + ...], (21) 

while in the opposite regime, when A S> y/Dot/A, the density past the rightmost particle approaches 
the initial value p exponentially fast. It is important to notice that Eq.(21) shows that the density past 
the rightmost particle is almost constant (and different from p) in a certain region whose size grows 



12 



in proportion to X(t). We note finally that for the just described one-dimensional model Eqs.(17) to 
(19) are exact, as it was shown subsequently by rigorous probabilistic analysis in [56]. 



Now, Eq.(18) predicts that four different regimes can take place, depending on the relation between 
p and p. First, for p < 1 — p the parameter A is negative and thus the rightmost particle effectively 
compresses the gas phase. When p exactly equals 1— p, (which defines the "yield" value of p necessary 
to initiate further compression), the prefactor A in Eq.(17) is exactly equal to zero, A = 0. Thus X(t) = 
and the gas phase is stable. In this regime, however, despite the fact that X(t) = 0, the rightmost 
particle still wanders randomly around the equilibrium position and its mean-square displacement 
X 2 {t) grows with time. In |j55| it was shown that the growth is sub-diffusive and X 2 (t) ~ (1 — p)t 1 / 2 / p. 
Further on, A is positive and finite when 1 < p < 1 — p holds; in this regime entropic effects overcome 
the pressure exerted by the rightmost particle and the gas slowly decompresses. Finally, when p = 1, 
the rhs of Eq.(18) diverges, which means that A is infinitely large in the steady-state. Actually, in 
this case A shows a slow, logarithmic growth with time. At sufficiently large times, 

A « ^2 In(^) (22) 

Turning now back to the system under study, we have to identify p and to express it through 
the strength of interactions Uq(T), initial coverage p and the temperature As we have already 

mentioned, the asymmetry in the hopping probabilities of the edge particles or, in other words, the 
"edge" tension force arises from the attraction of the edge particles to the "bulk" monolayer. Thus, 
since the interaction potential rapidly vanishes with the distance, we may expect that it is mostly 
dominated by the density profile in the vicinity of the moving edge, which is itself dependent on the 
magnitude of the edge tension. Self-consistent choice of p is thus prescribed by Eq.(14), which relates 
the hopping probabilities p and q to the density profile in the monolayer. Taking into account the 
result in Eq.(21), which defines the density profile in the vicinity of the moving edge, we find from 
Eq.(14) that the ratio p = p/q obeys 

p = exp(-Pa~/ edge ) (23) 
In Eq.(23) the parameter 7 e d 9e stands for the edge tension, 

le d9 e = (1 - P) (24) 

and 5 can be thought of as the number of broken cohesive bonds due to a hop away of the two- 
dimensional edge of the monolayer. Explicitly, 

where the shortenings f± stand for the vectors (X(t) ± <r, Y) and the sum extends over all lattice sites 
excluding r" = f±. For the simplest case of the square lattice, when z = 4, the parameter 8 can be 
readily calculated in explicit form: 

oo oo 

S = a-\2 + £ (1 + J 2 )' 3 } ~ 3.4 a~\ (26) 

j=l j=-oo 

which shows that in the presence of weak long-range attractive interactions 5 only slightly exceeds 
5 C = 3 a -1 - the result which we would obtain in the extreme case of nearest-neighbor (Ising-type) 
attractions. 
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Therefore, we have that the mean displacement of the monolayer edge obeys Eq.(17), i.e. X(t) = 
AyJ Dot, which is consistent with the behavior predicted in ^7| for spreading (A > 0) precursor 
films. In our case of a semi-infinite monolayer, the parameter A in this growth law is defined implicitly 
through Eqs.(18), (24) and (25), which allow to interpret it in terms of the strength of the PP 
interactions, initial mean coverage p and the temperature. In Fig. 4 we present numerical solution 
of these equations and plot A versus the dimensionless parameter e = f3Uo(T)a5/2, which appears to 
be the most significant critical parameter of the model. The solid curves in Fig. 4 define the function 
A{e) for four different initial densities. These curves demonstrate that the evolution of the monolayer 
is, in general, very sensitive to the value of e and it may proceed rather differently, depending on the 
relation between this parameter and p. 

We continue with some analytical estimates of the e-dependence of the prefactor A, where we can 
specify four different regimes. 

I. When e belongs to a finite interval < e < 1 (high temperatures or low PP attraction), Eqs.(22) 
and (23) possess only one trivial solution p = 1, which means that in this range of parameters the 
"edge" tension is exactly equal to zero and the monolayer thus spreads as a surface gas. The edge 
in this regime advances a bit faster that pure \/t-lsw and follows X(t) ~ (t ln{t)) 1 / 2 . The particle 
density past the edge is almost zero within an extended interval, which grows in proportion to X(t). 



II. In the range 1 < e < e c , where 

_ ln{l - p) 

P 



(27) 



the parameter A is finite and positive. Therefore, in this regime the monolayer also wets the substrate 
and the edge displaces in proportion to \/t. It is easy to check that in this regime the edge tension 
ledge is positive and vanishes as 

ledge ~ (T b - T) (28) 

when the temperature T approaches the value Tb. The critical temperature Tb is implicitly defined by 
the condition e = 1, which can be rewritten as 

n - M^l (29) 

Now, since the edge tension is positive below the T& and is exactly zero above the it seems natural 
to identify this regime as the regime of liquid-like spreading and, correspondingly, the temperature Tb 
- as the temperature of the surface gas-liquid transition or, in other words, as the boiling temperature 
of the monolayer on solid surface. We note finally that the parameter A diverges in the limit T — > Tb, 
(when e — > 1), 

A ~ sJh^) (30) 
Within the opposite limit e — > e c , the parameter A vanishes as 

A ~ (1 ~f )(£c ; e) (31) 
l-(l-p)e c 1 ; 

In this liquid-like spreading regime, the particle density past the edge is nearly constant within a 
region of size X(t) and is lower than the unperturbed density p in the bulk monolayer. 
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III. At the point e = e c the monolayer partially wets the substrate; p(\; t) = p and the prefactor 
A is exactly zero. We thus denote this point as the point of the wetting/dewetting transition. The 
corresponding critical temperature T w / dw is defined by Eq.(27), which gives, explicitly, 

_ Uo(T w / dw ) ° 8 . . 

J-w/dw — ^ ( 6Z ) 

Thus, this critical temperature appears to depend on the monolayer coverage p. We note now that 
two critical temperatures are simply related to each other. When Uq{T) is independent of T, like it is 
in the case of the London- van der Waals interactions, we find the following relation 

T w /dw = (33) 
For the Keesom-van der Waals interactions, when Uo(T) ~ 1/T, we find instead of Eq.(33), 

T w /dw = — f= 



For the monolayer, the wetting/dewetting transition point e c — > oo when p — ► 1; consequently, the 
critical temperature of the wetting/dewetting transition T w /d w — ► in this limit. Within the opposite 
limit, i.e. when p — * 0, T w / dw — > T&. 

IV. Finally, for e > e c , which corresponds to the limit of either low temperatures or strong particle- 
particle attractions, the parameter A is negative, A < 0, and thus the presence of a monolayer with a 
given coverage p on the solid surface is energetically non-favorable; consequently, it dewets from the 
substrate. The e-dependence of the parameter A appears to be very weak in this regime (see Fig. 4), 
which means that compressibility of the monolayer is very low, being strongly limited by the process 
of diffusive squeezing out of "voids" at progressively larger and larger scales. The density before the 
retracting edge is higher than the mean value p in an extended region which grows in proportion to 
X(t) (an analog of the rim in the hydrodynamic dewetting [33, 53, |35|, |36|), 



We hasten to remark, however, that the predicted weak e-dependence of the parameter A concerns 
only the situation in which thickening of the monolayer is not allowed and in which the particle motion 
can be viewed as an activated hopping motion between the wells created by the particle-solid interac- 
tions; particle-particle interactions are assumed to be small compared to the particle-solid interactions, 
such that they can be treated only as a small perturbation. For liquids in confined geometries, where 
the geometrical constraints themselves do not allow for the thickening of the monolayer, the process of 
squeezing of voids out of the "bulk" monolayer will be the only mechanism of the dewetting process. 
However, for sufficiently strong particle-particle attractions, comparable to the diffusive barrier Vi, our 
approximate description of particles dynamics is not justified; consideration of the edge tension as the 
only driving force will not be appropriate either. 

For monolayers on open solid substrates, thickening of the monolayer by forming progressively 
higher and higher "towers" |3^], represents an additional mechanism of the dewetting process, which 
may be, under certain conditions, more efficient than diffusive squeezing of voids. One can thus expect 
that for sufficiently strong particle-particle interactions the dewetting will be facilitated by thickening 
of the film, resulting in more pronounced e-dependence of the parameter A. As found in [32|, in 



this regime the mean displacement of the edge still follows the \/t-dependence, which means that A 
remains finite in the limit t —* oo. We may, however, only speculate about the e-dependence of the 
parameter A for such a process, since such a possibility is not included into the model. We sketch in 
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Fig. 4 a hypotetical behavior in this regime (the curve given by squares). Lastly, when Uq(T) becomes 
comparable to the adsorption barrier Ud, one may expect transition to the hydrodynamic dewetting, 
when the monolayer tends to form a macroscopically large droplet. This regime was examined (starting 
from sufficiently thick initial films, however) in |3^, |34|, |3^, |3(| and it is known that the edge in this 
regime displaces at a constant velocity, which means that here A does not tend to a constant value as 
t — > do, but rather increases indefinitely as time evolves, A ~ —t 1 / 2 . In Fig. 4 we mark the transition 
from the monolayer-dewetting regime to the hydrodynamic dewetting regime by the line of crosses. 
We can not, of course, identify precisely the value of the parameter e at which such a transition takes 
place; this calculation requires, again, elaboration of the model allowing for the thickening of the 
monolayer. 

We finally comment that our results interpret the notion of the two-dimensional volatility, which 
is commonly used in experimental literature on spreading of molecularly thin films, in terms of the 
parameters of the particle-particle interactions, particle density and the temperature. It is precisely 
the relation between the value of the parameter e, which is the measure of the particle-particle cohesive 
interactions, and the critical value e c , which shows whether having a monolayer on the solid substrate 
is energetically favorable or not. Consequently, we may expect that liquids with e > e c are not volatile 
in two dimensions, while liquids with e > e c are. 



6 Conclusions. 



To summarize, we have presented a microscopic dynamical description of the time evolution of a 
monolayer on solid surface. The monolayer was assumed to be created in initially non-equilibrium 
configuration, in which it covers only one half of the solid surface, and then was allowed to evolve 
in time due to particles random motion. Particles hopping motion was determined as the Kawasaki 
particle- void exchange dynamics in presence of long-range particle-particle attractions . Here we 
have focused exclusively on a two-dimensional behavior, assuming that particle evaporation from the 
substrate is absent and thickening of the monolayer is forbidden. We have shown that in such a 
situation the behavior of the monolayer is very sensitive to physical conditions and parameters of the 
particle-particle attractions. Depending on the strength of the latter, the monolayer can show different 
kinetic behavior; it can wet, partially wet or spontaneously dewet from the substrate. More precisely, 
our results can be summarized as follows. We find that the mean displacement of the monolayer edge 
X(t) evolves as X(t) = A(Dot) 1 / 2 , where Dq is the bare diffusion coefficient describing dynamics of 
an isolated particle on top of solid surface, while A is some parameter dependent on the strength 
of liquid-liquid attractions Uq{T) and temperature T. At sufficiently high temperatures, such that 
T > Xfc, the parameter A is greater than zero and shows a slow growth with time, A ~ {ln(t)) l l 2 . 
We identify this regime as "surface-gas" spreading, since the tension of the monolayer edge appears 
to be equal to zero. Actually, in this regime the density past the edge is almost zero in an extended 
interval, whose size grows in proportion to t 1 / 2 . Next, in the range of temperatures > T > T w i<i w 
the parameter A is also positive but tends to a certain constant value as t — > oo. The edge tension in 
this regime is positive and we thus call it as the regime of " liquid- like" spreading. The density past 
the edge is constant within the interval of size X(t) and is less than the density in the bulk monolayer. 
Thus in both regimes the monolayer expands and wets the substrate. Further on, at T = the 
parameter A = 0, i.e. there is no regular dependence of the displacement of the edge on time and the 
monolayer remains in its initial configuration. Finally, below the temperature of the wetting-dewetting 
transition, T w u w , the parameter A is constant and negative, i.e. the monolayer contracts by squeezing 
out the "voids" and dewets from the solid surface. 
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Fig.l. Initial configuration of a monolayer on top of solid surface. Wavy lines depict the 
potential energy landscape created by the solid atoms. 




Fig. 2. Liquids in confined geometries: (a) Uptake of liquid in a nanopore. (b) Ultrathin liq- 
uid film with a fluctuation-induced "hole" , sandwiched in a microscopically thin slit between 
two macroscopically large solid surfaces. 
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Fig. 3. Associated one-dimensional hard-core lattice gas model. Empty circles denote the 
hard-core gas particles, whose probabilities m of jumps to the right and to the left are 
symmetric. The filled circle stands for the rightmost particles with asymmetric hopping 
probabilities, p and q. 




Fig. 4. Numerical solution of Eqs.(18),(24) and (25). Solid lines from top to bottom show 
the dependence of the parameter A on e for p = 0.9, 0.8, 0.7 and 0.5 respectively. The inset 
displays behavior of e c (the point in which the curve A{e) crosses zero) as the function 
of the monolayer density p. Diamonds indicate the boundary line e = 1, which separates 
the" surface gas" and liquid-like behaviors. The crosses outline the demarkation line between 
the dewetting in the monolayer regime and the hydrodynamic dewetting (H.Dw), when 
thickening of the film can appear; the squares present the hypotetical continuation of the 
curve A(e) in the hydrodynamic dewetting case (see also explanations in the text). 
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